TNT 97 



Characterization of the spatial complex behavior and transition 

to chaos in flow systems 

M. Falcioni D. Vergni and A. Vulpiani WW 

Dipartimento di Fisica, Universitd "La Sapienza" , P.le A. Mow 2, 00185, Roma, Italy 
(*) Corresponding author (fax:+39-6-4463158; E-mail: M.Falcioni@romal.infn.it) 

(a) also INFN, Sezione di Roma 1 
(b) also INFM, Unita di Roma 1 
(February 5, 2008) 

Abstract 

We introduce a "spatial" Lyapunov exponent to characterize the complex 
behavior of non chaotic but convectively unstable flow systems. This com- 
plexity is of spatial type and is due to sensitivity to the boundary conditions. 
We show that there exists a relation between the spatial-complexity index we 
define and the comoving Lyapunov exponents. In these systems the transi- 
tion to chaos, i.e. the appearing of a positive Lyapunov exponent, can take 
place in two different ways. In the first one (from neither chaotic nor spa- 
tially complex behavior to chaos) one has the typical scenario; that is, as the 
system size grows up the spectrum of the Lyapunov exponents gives rise to 
a density. In the second one (when the chaos develops from a convectively 
unstable situation) one observes only a finite number of positive Lyapunov 
exponents. 
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I. INTRODUCTION 



The dynamical chaos is considered to be one of the main sources of complex behavior 
in a dynamical system. One of the main properties of dynamical chaos is the sensitive 
dependence of the evolution on the initial conditions, i.e.: a small error on the initial state 
grows exponentially in time This behavior is usually assumed as the characterizing 
property of chaos, and it is quantified by a positive value of the maximal Lyapunov exponent 
Ai. 

However, highly nontrivial behaviors can appear also in systems which are not chaotic 
(i.e., Ax < 0). Let us mention the systems with asymptotically stable fixed points, but with 
fractal boundaries of the attraction basins 0, and the chaotic scattering phenomenon j3|, 
where the "chaos" is just transient. 

An interesting situation can occur in high dimensional systems, like the following chain 
of maps with unidirectional coupling: 

x n (t + 1) = (1 - c)f a {x n {t)) + c/«(z n _i(t)), (1) 

where t is the discrete time, n = 1, 2, 3, . . . , N is a spatial index and xo(t) is a given boundary 
condition. These models are quite natural candidates for the description of flow systems, 
that are systems with a privileged direction, e.g.: boundary layer, thermal convection and 
wind-induced water waves [Q. 

After the seminal papers of Deissler and Kaneko || it is now well known that nontrivial 
phenomena can take place in systems with asymmetric couplings, even in the absence of 
chaos (Ai < 0). In particular, if the system is convectively unstable the spatial structure can 
be very complex and the external noise can have an important role in the formation and the 
maintenance of the structure [§,0. In spite of the clear evidence of a spatial "complexity" in 
these non chaotic systems, up to now, as far as we know, there is not a simple and systematic 
quantitative characterization of this phenomenon. To answer this purpose, we define a 
quantity that measures the degree of sensitivity of the system to the boundary conditions, 
and we study its possible quantitative relation with the comoving Lyapunov exponents - 
the quantities by means of which one can define the convective instabilities. The definition 
of the comoving Lyapunov exponents A(t>), for these extended systems, may be given as 
follows ||]. If 5xo(0) is a perturbation on the boundary at the time t = 0, in a frame of 
reference that moves along the system with velocity v > 0, at large t, this perturbation is 
O(Sx (0) exp[\(v)t}). If A(f) < for all v > the system is said to be absolutely stable; 
if there exists a range of velocity for which X(v) is positive, then the system is said to be 
convectively unstable. The interesting situation, in a convectively unstable system, arises 
when the usual Lyapunov exponent, Ai = A(t> = 0), is negative. 

In sect. II we discuss some results about the flow system model (JTJ) with f a (x) = 
ax (1 — x), that is the local map we use for all the computations. In particular we report 
on the qualitative spatio-temporal behaviors at varying the control parameters c and a. 

In sect. Ill we introduce an index - we call it "spatial Lyapunov exponent"- that 
supplies us with a quantitative characterization of the spatial complexity, in terms of the 
spatial sensitivity to the boundary conditions. We show that there exists a strong relation 
between the "spatial" Lyapunov exponent and the comoving Lyapunov exponents. 

In sect. IV the reader can find a study on the transition to chaos for the flow systems. 
The transition to a positive value of Ai can take place in two different settings: 
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I) from an absolutely stable state (i.e., a state for which the comoving Lyapunov exponents 

are negative for all velocities), that is not spatially complex; 

II) from a state that is convectively unstable (i.e., a state for which the comoving Lya- 

punov exponents are positive for some velocities), that has a certain degree of spatial 
complexity. 

The case I) is a rather standard transition, by which we mean that, for large N, there 
exists a finite density of positive Lyapunov exponents. On the other hand, in the case II), 
at varying N, one obtains only a finite number, not a finite fraction, of positive Lyapunov 
exponents. The latter behavior is a clear indication that the transition is dominated by a 
sort of boundary effect. 

Sect. V is devoted to conclusions and discussion. 

II. QUALITATIVE BEHAVIOR OF THE MODEL 

In this section we present some dynamical features of the unidirectionally coupled map 
lattice ([!]), with f a (x) — ax (1 — x). In many papers the boundary condition is kept fixed, 
i.e. xo(t) = x*, and often x* is an unstable fixed point of the single map x(t + 1) = f a (x(t)) 
0. Here, following Deissler [|| and Pikovsky ||, we adopt a more general time dependent 
boundary condition: x (t) = y(t) with y(t) a known function, that may be periodic, quasi- 
periodic or obtained by a chaotic system. 

At varying the control parameters, c and a, one observes a plethora of different spatio- 
temporal behaviors. In Fig. [I] we show the results of a numerical exploration of the phase 
space of the system, with a quasi-periodic boundary condition x (t) = 0.5 + 0.4sin(u;i), with 
to = 7r(v5 — 1). These results can be summarized by saying that, basically, there exist four 
qualitatively different spatio-temporal behaviors. 

A) Non chaotic and convectively stable. The comoving Lyapunov exponents are negative 

for all values of v (X(v ) < 0, Vt> ) and x n — > for n — > oo. The region of the parameter 
space corresponding to this behavior (absolute stability) is identified by the '□' symbols 
in Fig. P we call it "region A". One can say that the quasi-periodic (or chaotic) 
boundary condition x (t) is not able to excite the bulk of the system (see Fig. 0). 

B) Non chaotic and marginally convectively stable. The comoving Lyapunov exponents 

have a maximum value equal to zero for a v* ^ : \ max (v) = \(v*) = (the region 
with '+' symbols, in Fig. [1], th & t we call "region B"). The quasi-periodic boundary 
conditions produce spatio-temporal "strips patterns" (see Fig. |]). 

C) Non chaotic but convectively unstable. The maximal Lyapunov exponent, Ai = A(0), 

is negative, but the comoving Lyapunov exponent spectrum is positive in a certain 
interval of v > (the region with V symbols, in Fig. ffl, that we call "region C"). The 
spatio-temporal behavior appears irregular (see Fig. |j). 

D) Standard chaotic (the region with in Fig. |T], that we call "region D"). In this case 

one has A(0) > 0, and the spatio-temporal behavior is irregular, and is similar to the 
one of Fig. |j. 
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We remark that the results discussed above remain valid, with slight changes, for chaotic 
boundary conditions; e.g. with Xo(t) given by the y variable of the Henon map: y(t + 1) = 
—ay(t) 2 + (3y{t — 1) + 1, with typical values of the parameters a = 1.4 and (3 = 0.3. 



III. QUANTITATIVE CHARACTERIZATION OF THE SPATIAL BEHAVIOR 

We discuss here how to characterize the convectively unstable region C of Fig. |I[ Part of 
the results in this section have been briefly discussed in ref. 0. Some authors, e.g. Pikovsky 
H and Kozlov et al. fllPA , stressed the fact that the "irregularity" of these systems seems to 
increase with n. An analysis of function of t (by means of some standard methods 

for the characterization of dynamical systems, like, for instance, the power spectrum, the 
return map, the Grassberger-Procaccia correlation dimension |TT[) typically shows that x\ is 
more irregular than Xq, X2 more irregular than X\, and so on. Fig. |5| shows x n (t + 1) vs x n (t) 
for different n: it is evident an increasing of the irregularity as n increases. A simple way 
to characterize quantitatively the spatial complexity is by studying the spatial correlation 
functions: 

1 \ { x n x m} v^n) \*^m) /<->\ 

C{n,m) = — — — , 2 

where the average is with respect to the time. C(n, m) vs m, computed in the convectively 
unstable region, is shown in Fig. ^, for different n; two facts are evident: 

• the shape of C(n,m), at least for n ^> 1, does not depend on n, but only on \n — m\, 
thus revealing that we are observing a bulk property of the system; 

• the correlation decays exponentially, C(n,m) ~ exp(— \n — m|/£). 

It is natural to wonder how an uncertainty 5xo(t) = 0(e) - with e 1 - on the knowledge 
of the boundary conditions will affect the system. In this paper we consider only the case 
of infinitesimal perturbations, so that we may safely assume that Sx n evolves according to 
the tangent vector equations of the system ([!]): 

5x n (t + 1) = (1 - C)f a (x n (t))5x n (t) + cf'^Xn^tySXn^t). (3) 

For the moment we do not consider, for sake of simplicity, intermittency effects, that is, 
we neglect finite time fluctuations of the comoving Lyapunov exponents. The uncertainty 
Sx n (t), on the determination of the variable at the site n, is given by the superposition of 
the evolved 5x (t — r) with r = n/v: 

5x n {t) ~ j Sx Q (t - r)e x(v)T dv = e J e [x(v)/v]n dv. (4) 

Since we are interested in the asymptotic spatial behavior, i.e. the large n one, we can write: 

5x n (t) ~ ee^\ (5) 
where, in the particular case of a non intermittent system, one has: 
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7 = max . (6) 

v v 

Equation (|6]) is a link between the comoving Lyapunov exponent and the "spatial" Lyapunov 
exponent 7, a more precise and operative definition of which is given by: 

7=lim 1 6n^\ ( 7 ) 



n— >oo 



n \ e 



where the brackets mean a time average. 

So equation (|) establishes a relation between the convective instability of a system 
and its sensitivity to the boundary conditions, which can be considered a sort of spatial 
complexity. 

We remark again that equation @ holds exactly only in the absence of intermittency; 
in this case it can be shown from (|6]) that our spatial index can be written in a simple way 
in terms of another kind of spatial Lyapunov exponents, /i(A), introduced in ref. 



7 = max ^4 = //(A = 0). (8) 

In the general case the relation is rather more complicated. We introduce the effective 
comoving Lyapunov exponent, A t (u), that gives the exponential changing rate of a pertur- 
bation, in the frame of reference moving with velocity v, on a finite time interval t. 

Then, instead of (£|) we obtain 

5x n (t) ~ e J e [ ' Xt{v)/v]n dv, (9) 

and therefore: 

7=lim-(lnJ — -)=lim-lnJ — - -= max^ . (10) 

n^oo n \ e I rwoo n e \ v V I 



In a generic case, because of the fluctuations, it is not possible to write 7 in terms of X(v). 
Nevertheless it is possible to state a lower bound: 

7 > max = max = 7 . ( 11 ) 

The evaluation of the function A(t>) needs a heavy computational effort, however one can 
find good approximations of the quantity 7*. A first simple approximation for it, actually a 
lower bound, is given by 

n - (12) 

V* 

where v* is the velocity at which A attains its maximum value. The analysis of the long 
time behavior of many impulsive perturbations makes it possible to obtain v* and A(t>*) 
without the knowledge of A(v) as a function of v. An improvement of this approximation 
can be performed in the following way. Beside A(i>*), one computes the usual Lyapunov 
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exponent Ai = A(0), then one estimates the function X(v), by assuming it is the parabola 
\ p (v ) passing in the point (0, Ai) with maximum X(v*) for v = v* and, finally, one determines 
7 P = max v [X p (v)/v}. Typically 7 P is very close (within a few percent) to 7*. 

In Fig. [7] we show 7, 7*, 71 and 7 P versus a at a fixed value of c (c = 0.7), again for the 
logistic map with the quasi-periodic boundary condition xo(t) = 0.5 + 0.4sin((v^5 — l)%t). 
There is a large range of values of the parameter a for which 7 is rather far from 7*; for 
instance, at a = 3.74 we have 7 = 0.28 and 7* = 0.26. 

The difference is an effect of the intermittency; this may be pointed out by looking at 
what happens with the map f a (x) = ax mod 1: in this case we find that, all over the explored 
range of variation of a, 7 and 7*, from a numerical point of view, are indistinguishable (their 
relative difference is smaller than 10~ 6 ). 

We may obtain a further indication of the fact that the non negligible fluctuations of the 
comoving Lyapunov exponents are at the origin of the marked difference of 7 from its lower 
bound, by introducing, following ref. the generalized spatial Lyapunov exponents, L s (q). 
These quantities allow us to characterize the fluctuations in the growth of the perturbations 
along the chain: 



5x n 



(13) 



L s (q) = lim — ln( 

By means of standard arguments of probability theory, one has that: 

• L s (q)/q is a monotonic non decreasing function of q\ 

• dL s (q) I dq\ q=0 = 7; 

• L 8 (q) = 75 + \o 2 q 2 , for small q, where a 2 = lim n _ >0O ((ln|5x n /e| — ■yn) 2 )/n . 

The shape of L s (q)/q depends on the details of the dynamics, however L s (q) is fairly deter- 
mined - at least for small values of q - by the two parameters 7 and a 2 . The reason for 
having introduced this function is that one expects some relation between the fluctuations 
of the spatial-complexity index 7, and the fluctuations of the effective comoving Lyapunov 
exponents, and it is much more easy to compute the former than the latter. Fig. |8|, shows 
that, in the case of the logistic map, as we expected, the parameter a 2 (that is related to the 
variance of the spatial fluctuations) is small (large) in the region where 7* is a good (bad) 
approximation of 7. 

We stress that the definition (|7|) and the bound ( |TTD have a general validity. It is not 
difficult to understand that they are valid not only for 1-D flow systems, such as model ([[]), 
but also for continuous-time systems (as in the case of the asymmetric Ginsburg-Landau 
equation). As a matter of fact, all the arguments discussed above hold unaltered whenever 
one can introduce the comoving Lyapunov exponents. 

At the end of this section we want to note that there is not a simple relation between the 
correlation length £ and the exponent 7, such as, for instance, £ oc 7™ 1 . Indeed - in analogy 
with the case for the corresponding quantities (characteristic correlation time and maximal 
Lyapunov exponent) used to characterize the temporal behavior of the dynamical systems 
with few degrees of freedom - we do not expect a simple relation between £ and 7. 
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IV. TRANSITION TO CHAOS 



From Fig. [I] one can see that, for the system under investigation, there exist two routes 
for reaching the chaos: 

I) the way from non chaotic and convectively stable behavior, i.e. from region A to region 

D; 

II) the way from non chaotic, but convectively unstable behavior, i.e. from region C to 

region D. 

In this section we study with a certain detail this twofold way to chaos, by looking at the 
features of the Lyapunov exponents of the system. 

It is easy to understand that, for the system (|]) the computation of all the Lyapunov 
exponents is much easier than in the generic case. The origin of this lucky fact is in the 
triangular structure of the Jacobian matrix M[x(t)] that rules the evolution of the tangent 
vector 



5x(t+ 1) = M[x(t)]5x(t), 
where <5x = (Sxi, 6x2, ■ ■ ■ , Sxn), x = (x±, X2, ■ ■ ■ , xn) and 



MM 





V '■ 





(1 " c)f a {x 2 ) 

cfLM 






(1 - c)f' a (x 3 ) 









(14) 



/ 



(15) 



Since the product of triangular matrices is again a triangular matrix, all the Lyapunov expo- 
nents Ai, A 2 , . . . , Ajy (ordered, as usual, according to the decreasing values) can be computed 
in a simple way - without using the standard orthonormalization method of Benettin et al. 
T4| - from the quantities 



A, = ln(l-c) + Um±£H/aO*(n))l, 

n=l 



1,2, 



N. 



(16) 



The Lyapunov exponents {Xj} then are nothing but the {A,} after a reordering by decreasing 
values. We stress that, for the system under investigation, the computer time Tjy for the 
computation of all the Lyapunov exponents is proportional to the system size N: T N ~ iV; 
while in a generic map, with local coupling, one has T/v ~ iV 2 . 

The transition of type I) - that takes place in the zone of the region D close to the 
boundary with the region A - shows the behavior already observed in the maps with a 
generic symmetric local coupling [[L5[]. By looking at Fig. |9] one clearly sees that there exists 
a very well established thermodynamic limit for the Lyapunov spectrum. When iV — ► 00 
there exists a limiting function G(x) such that 



A, ~ G{i/N). 



(17) 



The existence of this limit entails various consequences. From (|T7j ) one can conclude that 
the number of the non negative Lyapunov exponents N is proportional to TV: 
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N ~ N. (18) 

So, by the Pesin formula Q, one obtains a finite Kolmogorov- Sinai entropy per degree of 
freedom, h: 

H 1 N r 1 

h=hm-= lim -E V(A») = / G(x)£(G(*))da;, (19) 

where if is the Kolmogorov-Sinai entropy and # is the step function. In addition, from eq. 
( PI) and the Kaplan- Yorke conjecture |l[] one infers that the information dimension dj of 
the attractor is proportional to N: 

di ~ N. (20) 

For the transition of type II) - taking place in the part of the region D close to the 
boundary with the region C - the features of the Lyapunov exponents are very different 
from those described above. Fig. shows that the behavior of the {Aj} do not follow eq. 
(PI); on the contrary one has 

Ai~F(i). (21) 

Therefore, in this case, the Kolmogorov-Sinai entropy and the information dimension of the 
attractor are not extensive quantities: H = 0(1) and dj = 0(1) ViV. Loosely speaking we 
can say that this transition can be described in terms of a finite layer. 



At the end of this section we want to stress the following point. In Fig. [10], as a matter 
of fact, we show Aj vs i since, in the case of the transition of type II), the ordering of Aj 
according to their decreasing values coincides with their ordering according to the label of 
the lattice site on which they are computed (see eq. ([l6|) ): A n > A n+1 (with the exception 
of very few sites close to the boundary); this means that the horizontal coordinate in Fig. 
|T0| is just the site label of the chain. Therefore one realizes that, in this case, the positive 
Lyapunov exponents are in correspondence with the sites close to the boundary: this gives 
further support to the idea that we are observing a kind of finite layer phenomenon. On the 
other hand, for the transition of type I), where one has a good thermodynamic limit, there 
is no correspondence between the sites on the lattice and the Lyapunov exponents. This is 
well evident from Fig. [□]. 



V. CONCLUSIONS AND DISCUSSION 



In this paper we characterized, in a quantitative way, the spatial complex behavior and 
the transition to chaos of flow systems. We have shown that in a non chaotic, but convectively 
unstable flow - where the convective instability induces a spatial sensitivity to the boundary 
conditions - it is possible to introduce an index (a sort of "spatial" Lyapunov exponent) 
for the quantitative characterization of this "spatial complexity" . Moreover, there exists a 
relation (a bound) between this spatial complexity and the comoving Lyapunov exponents. 

The transition to chaos can take place in two possible scenarios: either from a state 
that is both non chaotic and spatially non complex, or from a state that is non chaotic 
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but convectively unstable. In the first case one has a standard thermodynamic limit: the 
Kolmogorov- Sinai entropy and the information dimension of the attractor are proportional 
to the system size. In the second case one has a sort of layer phenomenon, where only a 
finite number of Lyapunov exponents - related to the sites near the boundary - are positive. 

We conclude noting that all the results above do not depend too much on the details 
of the used boundary conditions xo(t). Indeed we have that if xo(t) has a chaotic behavior 
(like, for instance, that of a variable of the Henon map) 7 and a 2 , as functions of a are not 
very different from the case with Xo(t) a quasi-periodic function. The same is true for the 
properties of the two kind of transitions to chaos. 
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FIGURES 



FIG. 1. The behavior of the system (|l|) in the space of the control parameters c and a: '□' = 
absolute stability; '+'= marginal convective stability; '*'= connective instability; '■' = absolute 
instability. 

FIG. 2. Evolution of a state of the system, with a = 2.94 and c = 0.7 (region A); the boundary 
condition is quasi-periodic: xo(t) = 0.5 + 0.4 sin (uit), with u> = ir(\^5 — 1). 

FIG. 3. Evolution of a state of the system, with a = 3.40 and c = 0.7 (region B); the boundary 
condition is quasi-periodic, as in Fig. [2|. For a better graphical effect only the configurations at 
even times have been reported. 

FIG. 4. Evolution of a state of the system, with a = 3.85 and c = 0.7 (region C); the boundary 
condition is quasi-periodic, as in Fig. ||. 

FIG. 5. Graph of [x n (t), x n (t + 1)] at different n; the values of the parameters are c = 0.7, 
a = 3.9 (region C); the boundary condition is quasi-periodic, as in Fig. 

FIG. 6. Spatial correlation function C(n, m) = [(x n x m ) — (x n )(x m }}/[(x' 2 l ) — (x function 
of m, for n = 100 (a), n = 200 (b) and n = 400 (c) with parameters c = 0.7, a = 3.85 (region C); 
the boundary condition is quasi-periodic, as in Fig. [2|. 

FIG. 7. 7 (+), 71 (x), 7 p (□) and 7* (O) y s a at fixed c = 0.7; the boundary condition is 
quasi-periodic, as in Fig. ^. 

FIG. 8. a 2 (+) and 7 — 7* (□) vs a at c = 0.7; the boundary condition is quasi-periodic, as in 
Fig. |. 

FIG. 9. Xi vs i for the system ([!]), with c = 0.06 and a = 3.58, in the cases ./V = 100 (dotted 
line), N = 200 (dashed line) and N = 400 (full line); the boundary condition is quasi-periodic, as 
in Fig. g. 

FIG. 10. Xi vs i for the system (0), with c = 0.27 and a = 3.70, in the cases N = 200 (full line) 
and N = 400 (dotted line); the boundary condition is quasi-periodic, as in Fig. ^. 

FIG. 11. Non ordered Lyapunov spectra Aj vs i for the system (||), with N = 200 c = 0.06 and 
a = 3.58, obtained with two different initial conditions (identified with 'o' or '★') 
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